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Abstract. Our MAMBO 1.2 mm blank field imaging survey of ~ 0.75 sqd has uncovered four unusually bright sources, with 
flux densities between 10 and 90 mjy, all located in the Abell 2125 field. The three brightest are flat spectrum radio sources 
with bright optical and X-ray counterparts. Their mm and radio flux densities are variable on timescales of months. Their X-ray 
luminosities classify them as quasars. The faintest of the four mm bright sources appears to be a bright, radio-quiet starburst 
at z ~ 3, similar to the sources seen at lower flux densities in the MAMBO and SCUBA surveys. It may also host a mildly 
obscured AGN of quasar-like X-ray luminosity. The three non-thermal mm sources imply an areal density of flat spectrum radio 
sources higher by at least 7 compared with that expected from an extrapolation of the lower frequency radio number counts. 

Key words. Galaxies: active - Galaxies: high redshift - Galaxies: starburst - Galaxies: individual: MMI1540+6605, 
MMJ1541+6630, MMJ1541+6622, MMJ1543+6621 - quasars: general - Submillimeter 



Oh 1. Introduction 



Blank field surveys with the SCUBA and MAMBO bolome- 
ter arrays were able to resolve a good fraction of the extra- 
galactic far-infrared (FIR) background radiation that was dis- 
covered by COBE (Pu getet alJ (l996) into point sources, re- 
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vealing a population of previously unknown dust-enshrouded, 
optically faint starburst galaxies at high redshi fts ( Smail et al. 
| 1997[ iHughes et alJ Il998l iBarger et ap Il998l iBertoldi et all 
BOOa iBlain et al.ll2002HChapman et alJl2005h . One aim of the 
MAMBO survey was to cover a large field, in order to bet- 
ter constrain the abundance and nature of the brightest mm 
background sources. By now, about 2700 arcmin 2 have been 
imaged to rms depths ranging from 0.5 to 3 mjy per 11 arc- 
sec beam jBertoldi et alJ2000HDannerbauer et afll2002L [2004 ; 
lEales et all2003HGreve et alJ2004 . Most of the MAMBO and 
SCUBA sources show faint radio (20 cm) counterparts. Their 
radio-to-mm spectral energy distributions (SEDs) are similar 
to starburst galaxies, which are dominated in the radio by syn- 
chrotron emission, and in the (sub)mm by a steeply rising ther- 
mal dust emission spectrum. 

The estimated infrared luminosities of the (sub)mm 
background sources imply star formation rates of order 
1000 M yr _1 , which to sustain would require them to be very 
massive systems, comparable to the massive spheroids in the 
local universe. Recent CO detections in seven SCUBA and 
MAMBO galaxies imp ly luminous and dy namical masses of 
order 10 11 M Q dNeri et all2003UGreve et alJ2005h . 

It has been a main motivation for the wide field MAMBO 
surveys to constrain the volume density of the most massive 
starburst galaxies at high redshifts. We here report first results 
of our shallow wide field surveys, focusing on the detection of 
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Fig. 1. The MAMBO signal to noise 250 GHz map centered on 
Abell 2125. The rms noise level is below 1 mJy per 11 arcsec 
beam in the inner contour (ca. 400 arcmin-), and below 3 mJy 
within the outer contour (ca. 1200 arcmin 2 ). 

the four brightest sources ever found in mm or submm blank 
field surveys. 

2. Observations 

Using the Max-Planck Millimeter Bolometer (MAMBO, 
iKrevsa et all ill 9981 l2002l) 1 1.2 mm (250 GHz) 37- and 117- 
element arrays at the IRAM 30 m telescope we have im- 
aged four fields: the Abell 2125 field covers an area of 
about 1600 arcmin 2 (Fig. [TJ, the COSMOS field covers 
~500 arcmin 2 , the Lockman Hole covers ~400 arcmin 2 and 
the NTT deep field covers ~200 arcmin 2 , adding up to 
2700 arcmin 2 imaged to an rms noise level ranging from 0.5 
to 3 mJy per 1 1 arcsec beam. More than 60 significant (> 4<x) 
sources have been detected, most of which show a faint (< 
100;uJy) 1.4 GHz radio counterpart (Bertoldi et al., in prepa- 
ration). 

We discovered four previously unknown sources with 
S i.2mm > 10 mJy in this MAMBO survey (Tables [T] and 0. 
They are all located in the Abell 2125 field, which is cen- 
tered on the diffu se, low d ensity galaxy cluster Abell 2125 
at z = 0.25 dOwen et alJ Il999). Their mm flux densities 
were monitored with MAMBO pointed observations, show- 
ing considerable variability on timescales of months for the 
three brighter sources MMJ1540+6605, MMJ1541+6622 and 
MMJ1543+6621 (Fig.|2j. Those were also observed simulta- 
neously at 1 mm and 3 mm with the Plateau de Bure interfer- 
ometer (PdBI), showing them to have spectra declining with 
frequency (Fig. 0}. 

De ep (6.5 yidy rms) VLA 20 cm images dOwen et alJ 
l2005al) show likely radio counterparts to all four sources, 
providing positions to about 0.5 arcsec accuracy. Short ob- 
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Fig. 2. 250 GHz flux density measurement of the four brightest 
sources in the MAMBO fields, plotted as a function of time. 
The measurements are slightly offset horizontally for clarity. 



servations were also made a 4.86 and 8.46 GHz in August 
200 with the VLA D configuration, which are reported by 
lOwen etafl J2005bl) . and 324 MHz observations were made 
in April 2002 (Owen & Clarke, private communication). 
We performed near-simultaneous cross-scan observations of 
MMJ 1540+6605, MMJ1541+6622 and MMJ1543+6621 with 
the Effelsberg 100 m telescope at 4.85, 8.35, 10.45 and 
14.6 GHz on March 24 th , 2005, and cross-scans at 1.41 GHz 
were taken at Effelsberg two month earlier. These near- 
simultaneous observations show fiat radio spectra. Further 
cross-scan measurements were taken at 4.85 and 10.45 GHz 
nine and eight month earlier, showing radio fluxes that are vari- 
able on timescales of months for all three sources. 

We found 1 .4 GHz flux measu rements f or the three brighter 
sources also in the NVSS catalog dCondon et alj 2002). fo r two 
of the m at 325 MHz in the WENSS catal og jRengelink et all 
I l997[) . and at 5 GHz in the 87GB catalog (Gregory & Condon 

All four sources are covered by our mapping of the 
Abell 2125 MAMBO field with the Spitzer Space Telescope, 
using the Multiband Imaging Photometer for Spitzer (MIPS) in 
photometry mode at 24 ^m (rms ~ 1 1 fj.]y) and the Infrared 
Array Camera (IRAC) at 3.6 fim (rms ~ 0.5 /Jy) and 5.8 /vm 
(rms ~ 3 yuJy). Three of the sources are also covered by 
the IRAC 4.5 and 8.0 yum (rms ~ 0.6 and 2.5 yuJy) paral- 
lel observations. The basic calibrated data have been further 
reduced using the mopex software package. A median filter 
was applied to substract the background from the images. The 
source flux densities were measured by fitting a point response 
function that was determined from bright sources in the field. 
MMJ1540+6605, MMJ1541+6622 and MMJ1543+6621 were 
detected in all bands, while MMJ 154 1+6630 was not detected 
in any band. 

Spectroscopic redshifts for M MJ 1540+660 5 and 
MMJ 154 1+6630 are not available. Mille r et alJ J2004 
report that MMJ1541+6622 has a redshift of 1.382 based on 
the C III] and Mg II lines. These broad lines indicate that the 
object is a quasar. 
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Fig. 3. The radio-to-IR SED of MMJ 154 1+6630. The dotted 
line shows a greybody with 7d us t = 50 K and (3-1.5 and 
a synchrotron spectrum with a = 0.7 at z = 3, whose radio- 
to mm flux ra tio is set acco rding to the correlation found in 
local galaxies JCarilli & Yu nil 19991) . The solid line is the SED 
of Arp 220 shifted to z — 3.5 and scaled to match the radio and 
mm measurements. 



A long slit spectrum of MMJ 1543 +6605 was obtained us- 
ing Gemini-N and GMOS on June 30 th , 2003 for a total expo- 
sure of 30 minutes. The B600 grating was used which produced 
a spectrum with about 3.3A resolution. The spectrum shows 
one strong wide emission line (Fig.|5} with a number of associ- 
ated pairs of absorption lines. This combination and the wave- 
length separation of the absorption doublets clearly indicates 
that the line is Mg II. The redshift we derive is 1.2362. Also 
the low level emission near 6600 A is very consistent with the 
Fe II emission line complex at the same redshift. 

We estimate a photometric redshift of for 
MMJ1541+6630 (Fig. [3}, using the correlation between 
the radio and FIR lumin osities found in the local s tar-forming 
galaxies dCondonll992UCarilli & Yunll 9991120001) . 

We found X-ray counterparts for MMJ1540+6605, 
MMJ 15 4 1+6622 and MMJ 1543 +6621 in the ROSAT 1WGA 
catalog (Whit e et all 1 1994 . MMJ1541+6630 is detected in a 
20 ksec ROSAT PSPC observation. MMJ 1540+6605 is de- 
tected in a 83 k sec archival Chandra observation (Fig. |6j 
IWang et all ( 12004 1. The 0.5-10 keV flux measured by Chandra 
is (1.1 ± 0.1) • 10~ 13 erg cirT 2 s -1 . 

The offsets of all identifications to the VLA positions are 
less than the l<x positional error quoted in the catalogs, ex- 
cept for the NVSS counterpart of MMJ1543+6621 and the 
ROSAT counterpart of MMJ 154 1+6630, for which the offsets 
are within the 2<x positional error. 



Fig. 4. The radio-to-IR SEDs of the three brightest MAMBO 
sources, arbitrarily scaled for clarity. The dashed lines connect 
the simultaneous PdBI measurements (90 and 230 GHz). The 
solid lines connect near-simultaneous Effelsberg measurements 
(below 20 GHz). The Effelsberg measurements were taken dur- 
ing the same day, except for the 1 .4 GHz measurements, which 
were taken two month earlier. For J 1540+6605 measurements 
below 8 GHz were not taken because of confusion with a neigh- 
bouring source. 
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Fig. 5. The GMOS/Gemini-N spectrum of MMJ1543+6605, 
showing the Mg II line at z = 1.2362. 



3. Discussion 

To quantify the significance of our source identifications, we 
calculated t he corrected Poissonian probability, P, of a ch ance 
association dBrowne & Cohenlll97 8: Dow nes et alJ ll986). For 
the identifications from the NVSS, WENSS and 87GB cata- 
logs, P < 0.1%, so that the found counterparts correspond 
to the MAMBO sources with > 99.9% confidence. For the 
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Table 1. The radio to mm source properties of the four brightest MAMBO sources 



Source (J2000) 


S 0.3 




S 1.4 


S 4.9 


S s 


S io 


S 15 


5 95 


S 250GHz 


z 


"95 




[mJy] 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 






154000.01+660551.6 


30.3 ± 


0.2 


24-26 


5-41 


22 ± 1 


21-28 


20 ±2 


26+1 


12-19 


unclear 


0.5 ± 0.3 


154137.19+663031.3 


<0.7 




0.10±0.02 


< 0.09 


< 0.054 








10±2 


Q + 1.9 

D -12 




154141.01+662237.9 


68-89 




39-89 


26-66 


28-62 


59-68 


52 ±2 


37±2 


11-28 


1.382 


1.4 ±0.3 


154321.31+662154.5 


50.3 ± 


0.4 


43-61 


72-121 


63 ± 1 


62-112 


56 ±1 


89±5 


29-92 


1.2362 


1.3 ±0.2 



Origin of the measurements: 0.3 GHz: VLA & Westerbork; 1.4 & 8 GHz: VLA, Effelsberg; 4.9 GHz: Effelsberg, VLA & Green Bank 300ft; 
10 & 15 GHz: Effelsberg; 95 GHz: PdBI; 250 GHz: PdBI & MAMBO 
The spectral index a is defined as S v oc v~ a in this paper. 

Table 2. The mid-infrared to X-ray properties of the four brightest MAMBO sources 



Source (J2000) 


S 24 




S 5.8 


S 4.5 


S 3.6/im 


R 


0. 1-2.4 keV 




Hardness 


A).l-2.4keV 










L"Jy] 


L"Jy] 


L"Jy] 


[mag] 


[erg cm -2 s" 


'] 




[erg s" 1 ] 




154000.01+660551.6 


1113 ±7 


498 ± 10 


380 ±8 


338 ±4 


301 ±3 


20.9 


(4 ± .6) • 10 


14 


1.0 ±.4 


(4.2 ± .6) • 


10 43 


154137.19+663031.3 


< 36 


< 7.5 


< 7.5 


< 1.2 


< 0.9 


> 27.2 


(3 ± .5) • IO- 


14 


0.4 ± .3 


(2.4 ± .4) • 


10 4S 


154141.01+662237.9 


450 ±5 


169 ±8 


103 ±5 


79 ±2 


53 ±2 


20.7 


CS ± .8) • 10- 


14 


0.6 ± .3 


(9.4 ± .9) • 


10 44 


154321.31+662154.5 


1371 ± 18 




217 ±6 




116 ± 2 


19.3 


(1 ±.1)- io- 


13 


0.4 ± .3 


(8.9 ± .7) • 


10 44 



Origin of the measurements: 24 jim - 3.6 jim: Spitzer MIPS & IRAC; R: KPNO; 0.1-2.4 keV: ROSAT PSPC 

Hardness: ROSAT HR1 is defined as (M + J - C)/(M + J + C) where M, J, C are the countrates in the 0.5-1.3, 1.5-2.1 and 0.1-0.3 keV bands. 
For the luminosities we assumed n M = 0.27, Q A = 0.73 and H = 71 km s _I Mpc" 1 . For MMJ1540+6605 we assumed a redshift of 0.5, which 
is suggested by optical photometry. 



ROSAT identifications the confidence is 93%, 83%, 98% and 
99% forMMJ1540+6605, MMJ1541+6630, MMJ1541+6622, 
and MM J 1543 +6621, respectively. 

As the absorption of X-rays by hydrogen decreases with 
increasing photon energy, the hardness ratios can be used for 
rough estimates of the absorbing column densities. Three of 
the sources show absorption by A^h ~ 1 — 4 ■ 10 20 crrT 2 , com- 
parable to the galactic value of 2.4 ■ 10 20 cirT 2 toward the 
Abel l 2125 field estimated from the Leiden/Dwingeloo survey 
llHartmann & Burtonl 199H) . For MMJ 1540+6605 we estimate 
A^h ~ (2 + 0.4) • 10 21 cm- 2 , which hints at intrinsic absorption. 

3.1. MMJ1 541 +6630: A thermal source 

MMJ 154 1+6630 has a faint 1.4 GHz radio counterpart with a 
flux density typical of the MAMBO or SCUBA galaxies. The 
1 .4 GHz flux density measurement and the upper limits at 4.9 
and 8 GHz indicate a falling radio spectrum, consistent with a 
typical synchrotron spectrum with a ~ 0.7 (S v oc v~ a ). The 
ratio between the 1.4 and 250 GHz flux densities is consis- 
tent with a typical starburst galaxy SED at z ~ 3. No opti- 
cal counterpart could be identified down to an R magnitude of 
27.2, which is not unusual f or even the brightest mm galax - 
ies, e.g. MMJ 120456-074 1.5 JPannerbauer et alJl2002l Eool . 
The standard interpretation of these data is that this object is a 
luminous high redshift starburst. 

In starforming galaxies in which no AGN contributes sig- 
nificantly to the X-ray emission the radio and the FIR luminos- 
ity a re also found to be well correlated with the X-ray luminos- 
ity (Eanalli 2004). The X-ray source close to the radio position 
of MMJ 154 1+6630 is detected at a low signal to noise ratio, 



making the position determination uncertain. The confidence 
of an association of the X-ray source and MMJ 154 1+6630 is 
83%. In the following we assume that the ROSAT source is 
associated with MMJ 154 1+6630. The inferred X-ray luminos- 
ity is ~ 80 times higher than that expected from the correla- 
tions between X-ray, radio and mm luminosities for starform- 
ing galaxies. This X-ray excess reveals the presence of a radio- 
quie t AGN, which is a lso not unusual for (sub)mm galaxies: 
Alexande r et alJ l|2004) suggest that ~ 40% of the (sub)mm 
galaxies host AGN, which however contribute less then 20% 
to the IR bolometric luminosity. Fol lowing the classification 
scheme suggested by Hasinger (2003), the implied X-ray lumi- 
nosity classifies MMJ 154 1+6630 as a quasar. 

Taking the SED of the local ULIRG Arp 220, redshifted 
to z = 3.5, and fitting it to the mm and radio measurements 
of MMJ 154 1+6630 would yield flux densities in the Spitzer 
IRAC bands that are 3.5 to 20 times higher than the measured 
IRAC upper limits. A typical ULIRG SED with emission from 
hotter dust (T > 50 K) than that in Arp 220, located at higher 
redshi fts than 3 wo uld still match the radio and mm measure- 
ments felainll 19991) but would yield lower flux densities in the 
IRAC bands. It could be that MMJ 154 1+6630 is therefore at 
z>3. 

Scaling the a verage SED for a radio-quiet quasar from 
lElvis et alTll 1994b to match the observed X-ray flux, we would 
expect MMJ1541+6630 to have an R magnitude of 20.5, much 
brighter than what is observed. The expected flux densities in 
the Spitzer IRAC and MIPS 24 pun bands would also exceed 
the measurements by an order of magnitude. Dust obscuration 
by > 6.7 magnitudes at A = 0.16 fim (rest frame) could ex- 
plain the optical faintness. Assuming galactic dust properties 
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this would translate to an Ay > 2.6. To match the flux densities 
observed in the Spitzer bands, with galactic dust properties Ay 
needs to be of order 50 - 100, which is not atypical for star- 
burst galaxies; e.g. in Arp 200 Ay ~ 50. Assuming a galactic 
gas to dust ratio of 150, Ay > 50 would imply a hydrogen col- 
umn density Nh > 9 ■ 10 22 cirT 2 , which is about two orders of 
magnitude larger than the value inferred from the X-ray hard- 
ness. Due to the low signal to noise ratio of the ROSAT mea- 
surement, the X-ray inferred column density is uncertain by at 
least a factor of five. Another uncertainty arises from the high 
redshift, whereby the soft part of the X-ray spectrum, which 
is most affected by absorption, is shifted out of the ROSAT 
bandpass, so that much higher column densities would still be 
consistent with the observed hardness ratio (see e.g. Alexander 
et al. 2004). The difference in inferred column densities might 
also hint on different line of sights toward the starburst and the 
quasar. 

Studies of quasars at mm wavelengths show almost no cor- 
relation b etween the optical e mission of quasars and their mm 
emission dOmont et alJl2003l and references therein). Radio 
quiet quasars studied at cm and (sub)mm wavelengths show the 
same ratio of cm and (sub)mm inferred luminosities as known 
from starforming galaxies (Petric et al. in preparation; Beelen 
et al. in preparation). Besides this indirect evidence for a sig- 
nificant, if not dominant dust heating by a starbur st in quasars, 
for th e high redshift quasars PSS J 2322+1944 dCarilii et all 
l2003h . SDSS Jl 14816.64+5 25150.3 dBertoldi et alJl2003l) and 
the Cloverleaf (H1413+1 17. lweifi et all2003l) a starburst could 
directly be identified as the dominant dust heating source. The 
Spitzer measurements of MMJ 154 1+6630 rule out a signifi- 
cant contribution to the dust heating by an AGN, except if this 
source were a heavily obscured z ~ 5 quasar with moderate star 
forming activity. 

With these uncertainties in mind, we find the properties of 
MMJ 154 1+6630 consistent with typical X-ray luminous (sub- 
)mm galaxies such as those studied by Alexander et al. (2004). 
In the following we shall assume that MMJ 154 1+6630 is a 
quasar with its FIR emission dominated by star formation. 

For a dust temperature of 50 K and an emissivity i ndex 
B — 1.5, typical for local starburst galaxies dDunne & Ealea 
1200 lh . at z = 3 MMJ 154 1+6630 would have an infrared lu- 
minosity of (4 ± 1) ■ 10 13 L , making it one of the most lu- 
minous starbursts known. It would be more luminous in the 
(sub)mm and X-rays than the type-1 quasar SMMJ0 4135 that 
was r ecently discovered in SCUBA blank fields iKnudsen et alJ 



200 31), or the tvpe -2 quasar SMMJ02399 dlvison etalJ ll998: 



Gen zel etalJ Eo03). Note that the relation between 250 GHz 
flux density and infrared luminosity is nearly independent of 
redshift from z ~ 0.5 to 10. 

With a Sal peter initial mass f unction and a starburst age 



of ~ 50 Myr (O mont et alJl200lb the inferred star formation 
rate would be 4500 M yr _1 . Such a high star formation rate 
would imply a very massi ve object. Ad opting a dust emissivity 
Ki25^m = 18.75 cm g 1 (Hildebrand 1983) the inferred dust 
mass of MMJ 154 1+6630 is - 1 • 10 9 M o . 

The likely detection of a massive starburst associated with a 
quasar adds to the growing evidence that the growth of central 
supermassive black holes and the formation of the first stel- 
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Fig. 6. Smoothed ROSAT PSPC map of the Abell 2125 field 
with inlayed smoothed Chandra ACIS-I map. The logarithmic 
grey scale ranges from 1 ■ 10~ 3 ct/s to 6 • 10~ 3 ct/s in the ROSAT 
image and from 1 • 10 ct/s to 8 • 10~ 3 ct/s in the Chandra inlay. 
The positions of the four bright MAMBO sources are indicated 
by squares. 



lar po pulations are closely linked (e.g. lHaehnelt & Kauffmannl 
2000, and references therein). Today we see the fossil records 
of this connection in the correlation b etween black hole and 
stellar bulge masses in local galaxies dMagorrian et al] [l998: 
i Ferrarese & Merritd 120001 iGebhardt et all l2000h . Following 
lOmont et al.N200ll) . for MMJ 154 1+6630 we estimated the ra- 
tio between star formation rate, Msf, and black hole growth 
rate, Mbh, as well as the bulge to black hole mass ratio, 
Mbui/Msu. The former is estimated by relating the FIR lumi- 
nosity to the bolometric luminosity arising from the quasar, 
which we calculate from the X-ra y luminosity usin g the av- 
erage radio-quiet quasar SED from El vis et alTd 19941) . For the 
latter estimate we assume the dust properties mentioned above. 

We find Msf/Mbh ~ 3400. This is ~ 10 times higher than 
expected if black holes and stars formed simultaneously with 
the same ratio of bulge mass to black hole mass as observed in 
local galaxies. For SMMJ02399 we find M SF /M BH ~ 2200, 
still a factor of ~ 6 higher. For M bn i/M BH we find ~ 1200 
and ~ 850 for MMJ1541+6630 and SMMJ02399 respectively, 
a factor of ~ 2-3 above the value found in local galaxies. 
Given the uncertainties in our Mt, u i/M BH estimate and the scat- 
ter in Mbui/MeH for local galaxies, the estimated ratios are sim- 
ilar to those in local galaxies. Compared with the optically 
selected quasars of lOmont et al.l d200ll) . the bolometric lumi- 
nosities of the (sub)mm selected quasars MMJ 154 1+6630 and 
SMMJ02399 appear to be dominated by the starburst activity. 
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3.2. Flat spectrum sources 

The flat and time-variable radio to mm SEDs of the three 
brightest sources, their relatively bright optical counterparts 
and their high X-ray luminosities indicate that these sources are 
flat spectrum, radio loud quasars. The broad emission lines in 
the optical spectrum of MMJ 154 1+6622 also indicate that this 
source is a quasar. The Spitzer SEDs of MMJ 154 1+6622 and 
MMJ 1543 +6621 can well be fit w ith the average SEDs of ra- 
dio loud quasars jElvis et afll 19941) shifted to the spectroscop- 
ically measured redshifts. The 95 to 230 GHz spectral index 
measured with the PdBI also suggests a non-thermal spectrum, 
which is however steeper by ~ 0.5-1 than the radio spectrum 
between 1 and 95 GHz. The fact that the radio spectrum is flat 
up to 95 GHz indicates that the steepening occurs very close 
to this frequency. A similar mm-bre ak is also seen e.g. in 3 
quasa rs selected from the 3C catalog Jvan Bemmel & Bert oldi 
l200ll) . 

Flat spectrum radio sources (FSS) are well studied at ra- 
dio wavelengths, but little is known about their SEDs at fre- 
quencies above 90 GHz and below the MIR. To see how the 
mm bright FSS relate t o the radio selec t ed FS S, we used the 
90 GHz observa tions of | Ho ldawav et alJ fl*9 94). as interpreted 
by the model of Holdawav & Owen (2005), to extrapolate the 
source counts up to 250 GHz. This model uses 90 GHz obser- 
vations of flat spectrum quasars selected at cm wavelengths to 
sample the spectral index distribution of the sources' compact 
cores from 8.4 GHz to 90 GHz. To push the model to frequen- 
cies above 90 GHz, Holdaway and Owen determined plausible 
distributions of turnover frequencies in these FSS, assuming 
the spectral index turns over to 0.8 above these frequencies. 

The sampling of the spectral index distribution of the FSS 
popula tion up to 90 GHz builds on the work of IPatnaik et alJ 
(1992J), who selected sources with S schz > 200 mJy and 
1.4 t o 5 GHz spectral index a < 0.5. The 1.4 GHz 
data JCondon & Brodericklll985l Il986l) and the 5 GHz data 
(Condon et al. 1989J) that defined their su rvey w ere taken with 
the NRAO 91m telescope. IPatnaik etall Jl992h observed their 
FSS sa mple with the VLA A array at 8.4 GHz. lHoldawav et alJ 
( 1994) defined a subsample of Patnaik's sample, selecting 367 
sourc es with 8.4 GHz core fluxes ranging from 200 mJy to 
35 .Tv. lHoldawavetalJ i 19941) observed these sources with the 
NRAO 12 m telescope. As the 8.4 GHz observations with 0.2 
arcsecond resolution provided excellent core fluxes, and from 
8.4 GHz to 90 GHz, the steep spectrum emission is decreased 
by a factor between 10 and 100, a good spectral index distri- 
bution for the FSS cores between 8.4 GHz and 90 GHz was 
obtained. No evidence was found for a spectral index distribu- 
tion depending on core flux. 

ICondonl (l984) presents source counts of FSS at 5 GHz, 
complete down to 10 mJy, from a variety of single dish ob- 
servations. These counts, which included contributions from 
both a flat spectrum core and steep spectrum extended emis- 
sion, were reduced by a mean core fraction of 0.8 and then 
scaled up to 90 GHz with the 8.4 GHz to 90 GHz spectral in- 
dex distribution. To scale the FSS counts to other frequencies, 
we need a model for how the spectral index varies at higher fre- 
quency. Such a model is hinted at by the distribution of spectral 
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Fig. 7. The cumulative source counts based on the MAMBO 
and SCUBA surveys. The rightmost three data points corre- 
spond to the FSS and were not included in the counts below 
10 mJy. The lines show the model fit (solid) of the thermal 
population's counts and its l<x-error (dashed). The dotted lines 
represent the extrapolat ion of the 5 GHz FSS counts using the 
model of Holdawav & Owen (2005). 



index between 8.4 GHz and 90 GHz, which can be decom- 
posed into two Gaussian components, indicating that 81% of 
the observed sources have "turned over" to a steep spectral in- 
dex by 90 GHz, and 19% of the observed sources remain flat at 
90 GHz. By combining the theoretical knowledge that flat spec- 
trum quasars have a spectral index of about 0.0, and that steep 
spectrum, optically thin synchrotron emission has a spectral in- 
dex of about 0.8, and by making reasonable assumpti ons on the 
spread about those spectral index values, iHoldawav & Ower] 
(2005) were able to solve for a distribution of turn over fre- 
quencies. This distribution of turn over frequencies can then be 
combined with the well-studied 5 GHz FSS counts and the core 
fraction correction to estimate the FSS counts at any frequency 
above 8.4 GHz. 

One thing remains undetermined in this model: what hap- 
pens to those 19% of still-flat spectrum sources above 90 GHz? 
An upper limit to source counts at 250 GHz can be obtained 
by assuming that these sources remain flat spectrum through 
250 GHz. A lower limit can be obtained by assuming they all 
turn over right at 90 GHz. A "best guess" can be obtained by 
continuing the best fit power law distribution of turn over fre- 
quencies up to 181 GHz, at which point all of the remaining 
flat spectrum sources have turned over. 

The extrapolated number counts are shown in Fig.0 with 
no spectral steepening up to 250 GHz and with the "best guess" 
model for spectral steepening. The FSS source counts from the 
MAMBO fields are a factor of 9^ 3 or 7+f higher than pre- 
dicted from the extrapolation with and without steepening, re- 



Hauke Voss et al.: Quasars in the MAMBO blank field survey 



7 



spectively. Lensing amplification is unlikely to be a cause for 
the observed source overdensity since the four sources are far 
from the inner 2x2 arcmin region of the galaxy cluster where 
mild lensing amplification could be expected. It could be that 
the number of FSS in the MAMBO field is high by chance. 
With Monte Carlo methods we simulated a field of 1000 square 
degrees with source counts that are consistent with the model 
and sources at random positions. We took 1600 independent 
small fields of the MAMBO survey size from the large sim- 
ulation and analysed the distribution of source counts for the 
small fields. We found that less than 0.4% of these artificial 
fields contain three or more sources brighter than 10 mJy. This 
discrepancy might also be due to a deviation from the found 
spectral index distribution for sources below the 90 GHz detec- 
tion limit of 75 mJy. Although the spectral index distribution 
was found not to depend on radio flux density, this is based on 
the brighter sources (Ssghz > 200 mJy) that were also detected 
at 90 GHz. It is possible that some fainter sources have flatter 
or more inverted spectra and thus contribute noticeably to the 
mm number counts. 

The four sources are among the first quasars to be discov- 
ered in mm or submm blank field surveys. Although based on 
small number statistics, they allow a first study of the contribu- 
tion of non-thermal sources in (sub)mm observations. Figure[7] 
shows the cumulative source counts for the MAMBO and 
SCUBA deep field surveys. We have fitted the counts with a 
model of a population of dust enshrouded starburst galaxies. 
The model is based on the Sa unders 60 //m luminosity function 
of local galaxies JSaunders etalJl990h which consists of a flat 
power law for low luminosities and a steep power law for lu- 
minosities above ~ 10 9 L . An additional exponential cutoff at 
~ 10 12 L , and a luminosity and density evolution with redshift 
were applied to match the (sub)mm number counts and their 
redshift distribution. The model was found to be consistent 
with the SED of the extragalactic FIR background. Including 
MMJ 154 1+6630 in the number counts, we find that they are 
in good agreement with the extrapolation of the fainter source 
counts. The three brightest sources however indicate a flatten- 
ing in the source counts at 250 GHz flux densities > 10 mJy. 
The mm background appears to change from a declining pop- 
ulation of high redshift starbursts to a shallower distribution of 
less distant, radio loud quasars. 
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